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ABSTRACT 

An enigmatic and rare type of young stellar object is the FU Orionis class. The members are 
interpreted as “outbursting,” that is, currently in a state of enhanced accretion by several orders of 
magnitude relative to the more modest disk-to-star accretion rates measured in typical T Tauri stars. 

They are key to our understanding of the history of stellar mass assembly and pre-main sequence 
evolution, as well as critical to consider in the chemical and physical evolution of the circumstellar 
environment - where planets form. A common supposition is that all T Tauri stars undergo repeated 
such outbursts, more frequently in their earlier evolutionary stages when the disks are more massive, 
so as to build up the requisite amount of stellar mass on the required time scale. However, the 
actual data supporting this traditional picture of episodically enhanced disk accretion are limited, 
and the observational properties of the known sample of FU Ori objects quite diverse. To improve 
our understanding of these rare objects, we outline the logic for meaningfully constraining the rate 


of FU Ori outbursts and present numbers to guide 
surveys. 

Subject headings: stars: pre-main sequence —stars: 

tion —circumstellar matter 

1. INTRODUCTION 

Pre-main-sequence stars are highly variable due to 
mechanisms operating in the stellar photosphere, in the 
magnetosphere and innermost regions of the circumstel¬ 
lar disk, or in the disk atmosphere. While the phe¬ 
nomenology has received attention since the work of Joy 
(1945) and was pioneered in the CCD era by Herbst 
et al. (1994), it is only in the last several years that 
large scale, quality, multi-year duration and moderate- 
cadence (e.g. Findeisen 2015; Parks et al. 2014, Rice et 
al. 2012) or multi-week and high-cadence (e.g. Cody et 
al. 2010, 2011, 2014) sampling of the photometric time 
series phase space has become commonplace. These re¬ 
cently available data sets have enabled more complete 
characterization of the diversity of young star variabil¬ 
ity than was possible earlier, and they serve to further 
our understanding of physical processes associated with 
young stars. 

A number of identified or hypothesized variability 
mechanisms are illustrated in Figure Both periodic 
and aperiodic phenomena are represented, with short 
(minutes and hours) to long (years and decades) duration 
processes causing a wide range of brightness changes. Re¬ 
cent monitoring programs have contributed to the quan¬ 
tification of variability types, amplitudes, and timescales. 
But depending on survey parameters such as cadence, 
duration, and photometric precision, only some of the 
plausibly operating physical phenomena are identihable 
in any given time series data set. Thus, the relative dis¬ 
tribution of young stars in the plane of amplitude vs time 
scale shown in Figure is unknown. 

1.1. Typical Young Star Variability Patterns 
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parameter choices in the analysis of time domain 
variables: T-Tauri, Herbig Ae/Be —stars: forma- 


At lower amplitudes and intermediate time scales, sig¬ 
nificant progress has been made in understanding ob¬ 
served young star variability phenomena. Simple pe¬ 
riodic modulation is attributed to rotation of photo- 
spheric inhomogeneities across the pre-main sequence 
stellar disk. For “cool starspots” the analogy is to similar 
phenomena on the much older Sun, but with signihcantly 
enhanced amplitudes (up to ^ 0.1 mag). In addition, 
there may be rotationally modulated “hot starspots” at¬ 
tributed to accretion stream footprints. Variability of 
a periodic nature has been well studied over the past 
several decades, mostly in connection with a desire to 
measure angular momentum evolution and its relation 
to disks and accretion. 

By contrast, aperiodic variability of a stochastic or 
perhaps only short time coherence nature, is less well- 
understood. This is in large part due to the lack of 
obvious correlation between specific physics and specific 
variability patterns, though aperiodic behavior is gener¬ 
ally associated with circumstellar rather than stellar phe¬ 
nomena. Disk presence, accretion, and outflow are impli¬ 
cated. Variability amplitudes can exceed 1 mag but the 
typical rms is about 0.1 mag (e.g. Grankin et al. 2007, 
2008 and Venuti et al. 2015 for optical, Carpenter et al. 
2001, 2002 for near-infrared; and Morales-Calderon 2011 
for mid-infrared study). Findeisen et al. (2013) identified 
moderate amplitude, intermediate time scale fading and 
brightening categories in optical data, with fades lasting 
weeks to years and likely caused by obscuring material 
in the few tenths to few AU range, and bursts lasting 
days to months and likely arising in disk processes that 
are also outside of the magnetospheric region. 

Recently, Cody et al. (2014), Stauffer et al. (2014, 
2015), and McGinnis et al. (2015) published high qual¬ 
ity lightcurves of young stellar objects with disks in the 
young cluster NGG 2264. The combination of preci- 
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Fig. 1. — The variability parameter space occupied by em¬ 
pirically recognized phenomena (solid curves) and theoreti¬ 
cally hypothesized mechanisms (dashed curves) that operate 
in young stars. Some phenomena produce relatively sym¬ 
metric flux variations (e.g. starspots), while others produce 
fading events (e.g. eclipses, variable extinction) or bursting 
events (accretion variations). This representation is intended 
to be schematic only, with any detailed use requiring quantih- 
cation of the definitions of “amplitude” and of “time scale”. 
The FU Ori stars occupy the upper right of the diagram. 

sion, cadence, continuous duration, and wavelength cov¬ 
erage, revealed exquisitely detailed photometric fluctua¬ 
tions that were used to identify the following broad cat¬ 
egories of variability: 

• Flux dips with periodic or quasi-periodic nature are 
interpreted as circumstellar dust passing through 
our line of sight as it orbits, probably near the in¬ 
ner disk. The dust is hypothesized in different sce¬ 
narios as either entrained in accretion or wind flows 
between the inner disk radius and the stellar photo¬ 
sphere (explaining “narrow dip” stars) or located in 
warps in the inner disk regions (explaining “broad 
dip” stars a.k.a. AA Tau analogs). Flux dips that 
are aperiodic but repeating also likely related to 
dust in the accretion flow. 

• Brightening events relative to a more steady (but 
still variable) flux level are intepreted as unsteady 
mass transfer from inner disks to stars. The accre¬ 
tion flow has either variable flux along, or variable 
penetration of, a roughly dipolar magnetosphere. 
These events are aperiodic and there is no corre¬ 
sponding periodic category. 

• Phenomena that are rooted in the stellar photo¬ 
sphere include periodic modulation due to temper¬ 
ature inhomogeneities (cool/hot spots) that rotate 
with the photosphere. Aperiodic events include 
discrete flares due to coronal-like magnetic activ¬ 
ity. 

The recent plethora of monitoring data is enabling 
more detailed and quantitative study of the light curve 
shapes and the variability time scales (e.g. Findeisen et 
al. 2015) and amplitudes. The overall direction of work 


in the field now, is to tie a mathematical description of 
typical variability patterns among young stars more di¬ 
rectly to the driving physics. 

1.2. Rare Burst and Rarer Outburst Behavior 

Young star bursts and outbursts are indicated in Fig¬ 
ure with amplitudes above ^1 mag and time scales 
longer than ~1 week. The statistics regarding these phe¬ 
nomena currently are quite poorly constrained. 

Young stars that undergo repeated outbursts with ~2- 
4 mag amplitudes and characteristic burst timescales of 
weeks to months are called EX Lup type^ variables (Her- 
big et al. 2001; Herbig 2008; Lorenzetti et al. 2009). 
V1647 Ori type outbursts are also repeating, with larger 
^4-5 mag amplitudes, and similar or perhaps longer 
characteristic timescales of up to a year or so (Aspin 
et al. 2009, 2006). Finally, the FU Ori stars undergo 
^4-6 magnitude brightness increases on time scales of 
a few months to years (Herbig, 1989) and then decay 
over at least decades to centuries (empirically) or mil- 
lenia (theoretically). Because of the long times scales, 
repeated bursts have not yet been observed among indi¬ 
vidual members of the FU Ori class. 

To explain the moderate-amplitude bursts, there are 
a number of flavors of instabilties and pulsational be¬ 
havior in circumstellar disks that have been suggested to 
lead to variable and possibly cyclic inward mass flow. 
Figure [l] also shows the relevant ranges in amplitude 
and time scale for various models. Magnetic reconnec¬ 
tion instabilities are modelled by Lovelace et al. (1995) 
and Romanova et al. (2002), Rayleigh-Taylor instabil¬ 
ities investigated by Romanova et al. (2004), diffusive 
instabilities proposed by Goodson & Winglee (1999) and 
Romanova et al. (2005), and viscous instabilities advo¬ 
cated by D’Angelo & Spruit (2010, 2012). The largest 
amplitude and longest time scale events are thought to 
be driven by inner disk instabilities and associated rapid 
mass accretion, as described in detail below. 

The FU Ori outbursts are the most extreme and sus¬ 
tained, while the other burst categories are less extreme 
but more frequent. Even the stochastically variable burst 
behavior mentioned above as “typical,” is also probably 
driven by nonsteady or unstable accretion from the inner 
disk to the star. Thus there is likely a broad continuum 
of burst behavior. Understanding the frequency of oc¬ 
curence of different types of accretion burst events is im¬ 
portant for improving our understanding of stellar mass 
assembly. However, the balance of influence among these 
different burst categories on pre-main sequence and disk 
evolution is currently unknown. For the lower amplitude 
and more frequent events, rates are currently unstud¬ 
ied. For the higher amplitude and less frequent FU Ori 
stars, the event rate has been estimated in various ways 
in previous literature as discussed below, though num¬ 
bers range over several orders of magnitude. 

We focus here on the infrequent, well-defined, large 
amplitude burst events exhibited by FU Ori stars. The 
difference between modelling these events and modelling 
the related EX Lup stars (shorter duty cycle and lower 

^ These objects are sometimes called “EXor type variables” but 
they are actually named after EX Lup. The EXor notation that 
has crept into the literature is unfortunate since there is in fact a 
star EX Ori, but it is a pulsating AGB star. 
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amplitude) and V1647 Ori type objects (shorter duty 
cycle but similar amplitude, with evidence of the con¬ 
tributions of both accretion enhancement and extinction 
reduction in their brightening) is that the FU Ori stars 
can be identified for decades after their bursts, whereas 
the other types of objects must be identified during the 
weeks to months to year that they are bursting. We leave 
such an effort to a possible future contribution. 

2. THE FU ORIONIS CLASS OF OBJECTS 
2.1. Phenomenology 

The FU Ori stars are a very rare group of young pre- 
main sequence objects. The prototype FU Ori and the 
two other defining members of the class, V1057 Cyg and 
V1515 Cyg have their outbursts summarized from an ob¬ 
servational perspective by Herbig (1977), while Larson 
(1980), Herbig (1989), Hartmann & Kenyon (1985, 1996) 
discuss more general issues. Vittone & Errico (2005) and 
Reipurth & Aspin (2010) present recent review material. 
FU Ori itself, the prototype of the class, was known as an 
irregular variable star associated with the A Ori region; 
it underwent a ~6 mag outburst in 1936 and has decayed 
only 0.5-1 mag (where the range indicates the uncertainty 
in the peak magnitude) since this time. V1057 Cyg was 
discovered as emission line object LkHa 190 by Herbig 
(1958) and was presumed a normal T Tauri star until it 
underwent a 5.5 mag outburst in 1969 thus becoming the 
second known FU Ori star; its decay is much more rapid, 
and it is now only 2 mag above quiescence. V1515 Cyg, 
also a known variable before its outburst, experienced a 
much slower rise taking a dozen years to reach peak flux 
around 1980, 4 mag above quiescence, and like FU Ori 
decayed only very slowly. In each of the above cases, 
with the source luminosity increasing by factors of hun¬ 
dreds to thousands, the FU Ori outburst illuminated a 
reflection nebula. 

In addition to their outbursts, the FU Ori objects are 
further distinguished among young stars, in that their 
optical spectra show unique signatures of early spec¬ 
tral type and low photospheric density, similar to F or 
G supergiants (Herbig 1977) while near-infrared spec¬ 
tra suggest M supergiants (Mould et al. 1978, Sato et 
al. 1992, Greene et al. 2008). The spectra are broad¬ 
ened by several tens to roughly one hundred kms“^, 
much faster than the typical absorption line widths of 
accreting T Tauri stars. FU Ori stars also show spectral 
features and profiles associated with inflow and outflow 
processes, such as classical P Gygni line profiles in Ha 
and purely blueshifted absorption in higher Balmer lines, 
as well as similar signatures in He, Na D, Ga H H&K, 
and KI, attributed to winds. X-ray emission is also de¬ 
tected (Skinner et al. 2006). The broad band SEDs of 
FU Ori stars (e.g. Gramajo et al. 2014) exhibit large in¬ 
frared excesses due to circumstellar material that is both 
thermally emitting and viscously dissipating energy from 
rapid accretion. Important for the disk instability inter¬ 
pretation is the finding by Gramajo et al. (2014) that 
80% of FU Ori disks have M^isk > O.IMq, more massive 
than any of the disks among the well-studied sample of T 
Tauri stars in Taurus, and that 90% have accretion rates 
dMdisk/dt > 10“® Mq year“^, higher than 95% of Tau¬ 
rus T Tauri stars, though for any individual source one of 
course must be cautious of reported model parameters. 

The current census of FU Ori objects numbers only 


two dozen at most, some still debatable. 

2.2. Physical Interpretation 

To explain the high-amplitude outbursting FU Ori 
stars, theoretical models invoking thermal, gravita¬ 
tional^, or magnetorotational disk instabilities, and com¬ 
binations thereof, are proposed (Bonnell & Bastien 1992; 
Bell & Lin 1994; Clarke & Syer 1996; Kley & Lin 1999; 
Armitage et al. 2001; Vorobyov & Basu 2005; Bo- 
ley et al. 2006; Zhu et al. 2009; Martin & Lubow 
2014). There are also suggestions of instability triggering 
through dense cluster or close binary or planet interac¬ 
tions (e.g. Pfalzner 2008; Reipurth & Aspin 2004; Lodato 
& Clarke 2004). 

The FU Ori broad-band SEDs are well-reproduced by 
the disk models of Zhu et al. (2007, 2008) in which an 
outburst is driven by gravitational instability and asso¬ 
ciated magnetorotational instability that extends out to 
several AU at an accretion rate of 10“® Mq year“^ (Zhu 
et al. 2010). For the thermal instability mechanism (Bell 
& Lin 1994), by contrast, the outbursting zone is con¬ 
fined to <0.1 AU at 10“® Mq year“^ but could reach 
1 AU if 1O“^M0 year“^; the SED fits of this model are 
less good (Bell et al. 1995), although the light curves 
are well matched. The enhanced mid-plane disk accre¬ 
tion rates can be compared to those typically inferred 
for classical T Tauri stars which are more like 10“® Mq 
year“^ (Gullbring et al. 1998). 

In the standard picture (e.g. Hartmann & Kenyon 
1996), accretion outbursts are expected to decline in fre¬ 
quency and perhaps amplitude with increasing source 
age. More recent studies of protostellar accretion history 
(Vorobov & Basu 2010; Stamatellos et al. 2012, Dun¬ 
ham et al 2014; Bae et al. 2014) include disk instability 
physics in various ways in simulations, which results in 
the promotion of fragmentation. Accretion burst behav¬ 
ior can be produced at ages typically well less than ~1 
Myr, in the form of short-lived 1.5 to 2 orders of magni¬ 
tude enhancement in disk accretion rates. Padoan et al. 
(2014), conversely, advocate for the importance of vari¬ 
able infall (rather than disk accretion) rates, which are 
prolonged in their turbulent cloud fragmentation model 
relative to the more isolated core / envelope / disk mod¬ 
els that rely on disk instabilities to build up the stellar 
mass. 

It is widely assumed that FU Ori stars are typical 
solar-type young stellar objects with central star masses 
~0.3-1.5 Mq and ages <1-2 Myr. The predominant view 
is that they are Glass I type sources, with a significant 
amount of the luminosity in their quiescent stages coming 
from envelope infall. Some FU Ori stars, however, are 
demonstrated Glass H type sources such as PTF lOqpf 
(Miller et al. 2011), which was a routine M3e classi¬ 
cal T Tauri star with little or no evidence based on the 
SED for envelope infall Empirically, therefore, the 

^ Though not widely recognized, the disk instability model for 
FU Ori stars was first proposed by Paczynski in 1975, as reported 
by Paczynski, 1978 and Trimble, 1976. 

® The Gramajo et al. (2014) study identifies 1/3 of known 
FU Ori stars as Class II, but detailed examination shows likely 
errors in both directions. Most of their Class II designations 
(BBW 76 excepted) have flat or rising SEDs, as expected for Class 
I and so-called flat-spectrum sources, and some of their Class I 
designations (e.g. HBC 722 = PTF lOqpf) include far infrared 
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episodic accretion behavior appears to continue into the 
more optically revealed, envelope-free phase of star for¬ 
mation and pre-main sequence contraction, so possibly 
beyond 1 Myr of age. 

3. ESTIMATING OUTBURST RATES 

The frequency of an astrophysical phenomenon is of¬ 
ten derived by considering: the number of detections of a 
particular category of object, the biases for or against de¬ 
tection as well as the total number of sources in a parent 
sample from which the categories could be established, 
information about a relevant time scale. With these val¬ 
ues in hand one can arrive at a rate of occurence. A 
well-understood numerator and a similarly well-defined 
denominator are required for rigorous statistical assess¬ 
ment. 

3.1. Previous Approaches to FU Ori Outburst Rates 

Herbig (1977) advocated from the relative number 
of FU Ori stars and regular T Tauri stars known 
at that time, that the FU Ori outburst rate was 
10“"^ year“^ star“^, and that the phenomenon was re¬ 
current. 

Hartmann & Kenyon (1996) estimated the FU Ori out¬ 
burst frequency assuming the Miller & Scalo (1979) star 
formation rate within 1 kpc of the Sun, i.e. 0.01-0.02 Mq 
per year, and also seem to assume an average mass close 
to solar, giving ~0.01 new low mass stars per year. For 
the census at the time of 5 known outbursts (FU Ori, 
V1057 Cyg, V1515 Cyg, V1735 Cyg, and V346 Nor) in 
60 years within 1 kpc, an integrated outburst rate of 0.08 
per year is calculated. Dividing one rate by the other, 
a typical low mass star thus has ~10 outbursts. Divid¬ 
ing the T Tauri disk lifetime of 1 — 10 x 10® years by 10 
outbursts gives one outburst every 10®“® years. 

A related McKee & Offner (2010) estimation assumes: 
the Fuchs et al. (2009) star formation rate within 1 kpc 
of the Sun, 8 x 10“® Mq per year, and an average mass 
of 0.5 Mq, to get 0.016 low mass stars per year capa¬ 
ble of outburst. These authors also suggest that 25% 
of the mass of a star is accreted through such enhanced 
accretion episodes. 

An updated census of outbursting objects by Greene 
et al. (2008) using the catalog of Abraham et al. (2004) 
tabulates 18 spectroscopically identified FU Ori stars 
within 1 kpc; Reipurth & Aspin (2010) provide a more 
recent compilation. Dividing the number of presently 
known FU Ori objects by the star formation rate above, 
a typical low mass star spends ~1250 years of its life in 
an FU Ori phase. Although consistent with the Hart¬ 
mann & Kenyon figure for a 1-century outburst length, 
the coincidence probably isn’t meaningful since the in¬ 
crease in the FU Ori census in the last 15 years is offset 
by the higher star formation rate assumed. 

The above statements are equivalent to saying that 
the rate of FU Ori outbursts in active star forming re¬ 
gions is approximately 10“® year“^ star“^, independent 
of whether all T Tauri stars go through such a phase or 
whether only a subset is prone to outbursts. However, 
since the population of FU Ori stars is incomplete, this 
rate is likely an underestimate. 

photometry that is significantly confused by other sources in the 
beam. 


Another line of reasoning is that a typical T Tauri 
star must have a lifetime average accretion rate of ~ 
10“® Mq year“^ to account for its estimated mass and 
age. This high average accretion rate combined with the 
consequently lower-than-expected luminosity of “proto¬ 
stars” led Kenyon et al. (1990) to argue that FU Ori 
stars were an important mechanism by which much of 
the time could be spent accumulating mass at low rates 
and relatively little time at higher rates. If, for exam¬ 
ple, the star alternates between a quiescent phase of ac¬ 
creting 10“^ Mq year“^ (higher than typical) and an 
outburst phase of accreting 10 “^Mq year“^ (also higher 
than what is thought typical), then one can account 
for all the star’s accretion if it spends 1% of its life¬ 
time in outburst. Given an assumed typical outburst 
lifetime of ~ 100 years, this implies an outburst rate 
closer to 10“^ year“^ star“^ is needed, or an even higher 
10“® year“^ star“^ if the more typical quiescent and 
outburst accretion rate values of 10“® Mq year“^ and 
~ 10“® Mq year“^, respectively, are adopted. Of note 
is that BBW 76 has been at maximum light since at 
least 1900 (Reipurth et al. 2002), and V883 Ori since 
at least 1888 (Gederblad 1946); similar to FU Ori, these 
stars are not rapidly fading. Thus the actual duration of 
FU Ori outbursts may be a factor of several longer than 
what is typically assumed. Under the logic above, the 
outburst rate needed to reconcile the T Tauri accretion 
rate vs age problem using the episodic FU Ori outburst 
model is inversely proportional to the assumed outburst 
lifetime. Longer duration bursts mean lower rates and, 
conversely, if we allow for some fraction of shorter dura¬ 
tion FU Ori outbursts with lifetimes of only ~ 30 years, 
then the needed outburst rate is higher by factors of sev¬ 
eral than the numbers above. 

Bae et al. (2014) have recently produced a distri¬ 
bution of accretion rates during the infall and early 
disk accretion phases up to 1 Myr, based on 2D hy¬ 
drodynamic models that include heating and cooling. 
Roughly 26% of the time spent at the infall rate of 
~ 10-5-75 year“^, 64% of the time spent at the qui¬ 
escent disk rate of ~ 10“^ ’^® Mq year“^, and 10% of 
the time spent in outburst at ~ 10“® Mq year“^ (their 
Figure 16). This results in an average accretion rate of 
~ 1.5 X 10“® Mq year“^ over the first <1 Myr, more 
than enough to accrete the stellar mass using the simple 
logic above. 

The burst rate in these new Bae et al. (2014) mod¬ 
els is significantly higher than in previous gravitational 
instability scenarios where, inferring from their Figure 
3, about 24 outbursts in the first 0.3 Myr and only 2 
in the next 0.7 Myr are suggested from the more stan¬ 
dard models, whereas from their Figure 10 the number 
of oubursts is higher by at least a factor of two. Gonvert- 
ing to an estimated FU Ori rate results in values around 
8 X 10“^ year“^ star“^ in the Glass 1 stage, when the en¬ 
velope dominates the disk extending to about 0.25 Myr 
by which point the disk mass dominates, then dropping 
to about 3 X 10“® year“^ star“^ into the Glass H stage. 

Considering the results of the numerous lines of rea¬ 
soning above, it is not an exaggeration to say that there 
is currently a factor of at least 30 range among existing 
estimates of the FU Ori outburst rate. 

3.2. How to Measure the FU Ori Rate Empirically 



5 



Hypothetical Rate Tq star ' year ') 


We assume that, in practice, the estimated outburst 
rate r, in outbursts per star per unit time, will be in¬ 
ferred from the observations Nobs, for example using 
r = Nobs/e Nyso We assume that the uncertainty 
associated with the measurement r will be evaluated us¬ 
ing the confidence interval formalism. We then ask how 
large a survey must be to have a high probability of de¬ 
tecting enough outbursts (given fluctuations around the 
mean true outburst rate rg) so as to produce a tight con¬ 
fidence interval around the measurement r. 

Given an observed number of events Nobs and an in¬ 
ferred outburst rate r, an observer can define a confidence 
interval r'^PP®''), where the lower and upper lim¬ 

its on the outburst rate r are set by the desired confidence 
level C and by Nobs ■ For a symmetric confidence interval 
these limits satisfy the following equalities: 

^(^obs < ^Obs) ^ BiN^^^,Nyso,r^PP^^At) = 
^(^obs < ^obs) ^ 


Fig. 2.— The survey size needed to provide a factor of two 
constraint on the FU Ori rate at 90% confidence (C=0.9), 
with a 90% chance of detecting enough events to establish 
the constraint (R = 0.9). The minimum sample size is plot¬ 
ted as a function of the hypothetical “true” outburst rate rg 
for surveys having time baselines from 1 to 100 years. One 
may reduce the needed survey size by choosing a longer time 
baseline, by admitting higher uncertainty than a factor of 
two, or by requiring a lower confidence than 90%. 

One can consider it a reasonable goal to determine 
the frequency of FU Ori outbursts to within a factor 
of two. In the calculation that follows, we evaluate an 
experimental design intended to constrain the rate to a 
factor of two or better at conhdence C, with a chance 
of success (i.e., a measure of the risk in carrying out 
the experiment) R. Both C and R have values in the 
range [0,1]. We assume a region of sky containing a 
number Nyso of young stars with disks (massive disks 
are required in the standard instability scenario for the 


where in general ^(nQj^g, n,p) is the cumulative distribu¬ 
tion function for a binomial distribution with n trials of 
probability p that realizes riobs events. B is monotonic 
with respect to all three arguments: holding the other 
two arguments fixed, B increases as increases, de¬ 
creases as n increases, and decreases as p increases. 

For the experiment to have provided a factor of two 
constraint, a symmetric confidence interval around the 
nominal measurement r lies entirely inside the inter¬ 
val [r/2,2r]. Thus we need to consider the case where 
rlower > ir and r'^PP®^ < 2r. Substituting into the two 
equations above and expressing r in terms of Nobs, we 
want to design the experiment so that it satisfies, with 
probability R, the inequalities 


outbursts) that is monitored over duration At. We fur- j- ^ i i ^ • • i. i j. 

., 4 .U j. iv j- j- i i he adiustable parameters m this experimental setup are 

ther assume: that outbursts occur at some average rate ,, i m i ^ i aj. 

in-a -1 i -1 i .1 i. j-u i i i the sample size Ay CO and the monitoring interval At. 

rg = 10 year ^ star that outbursts do not repeat ^ at 

on observable timescales (i.e. less than a century, essen- ® ^ 1 -^obs’ Ayso, ) is a decreasing 

tially a constraint on rgAt ^ 1 star“^), and that any function of N^^^, while B^N^^g, Nyso, | sNyso ^ 

outburst occurring in the monitored interval is detected. increasing function. Therefore, for reasonable values^ of 

The number of observed FU Ori outbursts is then well jVyso and C, both these criteria will be satisfied for suf- 

modeled by a binomial distribution ficiently high A^j^g. For a factor of two constraint at 

/AT\ C = 0.90, we typically need A Uo ~ 8 for a broad range 

P(fc|A,p)= /(l-p)^-'= of Nyso. 

b / All that remains is to find the probability that A^j^g 

which, in our case with probability p = rgAt, becomes will be high enough, as a function of Nyso- Since N^^^ is 

, . drawn from a binomial distribution with Nyso trials and 

PiNovent\Nyso,roAt) = ( ] (rgAt)^'’''"*(1-rgAt)^^®?)rdbability rgAt {not the observed rate r that appears in 

\Nevent/ the preceding paragraphs), this can be done by choosing 


The number of FU Ori events that are expected to oc- t’rial values of Nyso, finding the minimum number of 

cur is Nevent = Nyso rg At and the number detected ^obs that would give a factor of two constraint on r, 

would be Nobs = tNevent, where we assume e = 1 in ^nd then checking whether the probability of actually 


what follows. The probability of Uobs or more detections 4 „ , . v, ^ at n rn, 

\ • 1,1 • , 1 r an example oi unreasonable values, let JMyso = 

PiNobs > nobs) IS given by the appropriate number of condition can be satisfied only if 7V„bs = b but B(l, 2, i) = 

binomial terms, F’(fc| Ayso, fgAt). 0.9375, SO the second condition cannot be satisfied if C > 0.875. 




6 


observing events, given the true rate ro, is at least 

our risk tolerance R. In other words, we predict the 
results of a future statistical analysis, while integrating 
over all possible experimental outcomes to which said 
analysis would be applied. The expressions above are 
solved simultaneously for the two parameters A^Qi^g and 
A^YSO ■ smallest value of IVygo which a solution 
exists is the minimum survey size. 

The trade-offs in parameter space between event rate 
r, survey baseline At, and minimum survey size Nyso 
are illustrated in Figure The survey size required to 
probe a given outburst rate is approximately inversely 
proportional to the available baseline (i.e., the expected 
event count is approximately constant). A sample of ~1 
million young stars with disks is needed in order to probe 
outburst rates ~ 10“® year“^ star“^ within a 1 year sur¬ 
vey. For longer time baselines, the survey samples can be 
smaller and still probe the same outburst rate - only 10^ 
stars are needed for a decade survey and 10"^ for a cen¬ 
tury baseline in order to probe the same rate as above. 
We note right away that typical studies today report in¬ 
formation on, at best, several thousand accreting young 
stars, and so can not yet even begin to meaningfully con¬ 
strain the FU Ori rate on empirical grounds. 

Figure demonstrates the effect of considering differ¬ 
ent time baselines and outburst rates for samples un¬ 
der 10® stars, a reasonable upper limit to the census of 
“known” young stars within a few kpc. Samples larger 
than this are required if we aim to measure in a short sur¬ 
vey the true FU Ori outburst rate and it is much lower 
than 10“^ year“^ star“^. 

Figure [^illustrates the impact of altering some of the 
fixed parameters in the simulations above: the factor of 
two sensitivity to deviations from the particular outburst 
rate, and the required probability of a successful experi¬ 
ment. As expected, a larger survey, having smaller sta¬ 
tistical fluctuations, provides either a higher probability 
of meeting specihc precision goals for the rate, or better 
precision on the rate. The bottom panel shows that rel¬ 
ative to the factor of two constraint we have specihed as 
adequate, improving our knowledge of the FU Ori rate 
to a factor of 1.5 requires four times as much survey data 
while knowledge to only a factor of 3 requires three times 
less survey data. The top panel illustrates the relative 
insensitivity of the survey size requirement to the risk 
of obtaining a less precise result on the rate than the 
experimental design. In other words, there is little to 
be gained by increasing the chance that a given survey 
size will not detect enough FU Ori outbursts to obtain a 
reliable constraint on the FU Ori rate. 

However, even if one sacrifices on quality of knowledge 
(Figure 1^, it is impossible to make useful statements 
about the FU Ori rate with a survey of fewer than sev¬ 
eral thousand stars over fewer than several years (Fig¬ 
ures ^ . Considering larger samples, for a survey of 

10,000 stars, an outburst rate of ^ 10“® year“^ star“^ 
can be measured or excluded after a few years, which is 
a currently possible exercise, while 10“^’® year“^ star“^ 
can be probed after 100 years, also feasible now by in¬ 
cluding particularly deep historical photographic plate 
data. A total sample size above 10® stars would system¬ 
atically shift the lines to more rapid returns, as would 
relaxing the requirements of 90% confidence or factor of 




Time Baseline (years) 


Fig. 3. — Cuts across Figurethat illustrate relationships 
among parameters for survey sizes consisting of less than 10® 
young stars with disks. Top panel shows the outburst rate 
ro for various values of time baseline At, while bottom panel 
shows the time baseline At for various values of the outburst 
rate ro probed. Both panels assume the same 90% chance of 
constraining the outburst rate to a factor of two. 

two knowledge of the rate. 

Finally, in Figure we explore the dependence of our 
ability to know the FU Ori rate on what the event rate 
actually is. As expected, the rarer the occurence of FU 
Ori events, the looser the rate constraint from any given 
survey. The tip of the plot at low rates and poor degree 
of knowledge of the rate, corresponds to a 90% chance 
of detecting at least one FU Ori event in the survey (in 
this case 10^ stars sampled over 30 years). 

4. APPLICATION TO MODERN SURVEY DATA 
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Experiment Success Probability Given Nqbs 



Factor to which Rate is Measured 

Fig. 4.— For a 30 year baseline and an outburst rate of 
10“^ year“^ star“^ (i.e. roAt = 3 x 10“^), the survey size 
needed to detect a factor of two deviation from the hypothet¬ 
ical outburst rate, at 90% confidence (O = 0.9), and 90% 
of the time {R = 0.9). Required survey size is plotted as 
a function of (top panel) the probability R that the survey 
will achieve the confidence bound, and (bottom panel) mea¬ 
surement precision on the hypothetical outburst rate. Larger 
surveys provide either higher probability of meeting specific 
precision goals on the rate, or better precision on the rate. 
Approximately 3500 stars are needed for factor of two preci¬ 
sion on the rate (bottom panel) at 90% confidence (top panel) 
for this time baseline and true outburst rate. 

Increasing interest in young star variability mecha¬ 
nisms has led to more and larger (both wider field and 
longer duration) as well as often multiwavelength time 
series data sets. The value of the earliest possible first 
epoch observation is considerable, as illustrated in Fig- 



5.0x10“=^ 1.0x10“^ 1.5x10-* 2.0X10“* 

True Event Rote {§ star“' year"') 

Fig. 5. — Dependence of our degree of knowledge of the 
FU rate on the rate itself. Simulation results are shown for 
a survey of 10"* stars over a 30 year baseline (same as Fig¬ 
ure assuming an experimental goal of obtaining a factor 
of two constraint on the FU Ori rate at C =90% confidence, 
and that we tolerate an R =90% chance of detecting enough 
sources to actually get the required constraint (as elsewhere 
in this paper). At tip of the plot where the constraint is 
reduced from a factor of two to nearly a factor of 5, the prob¬ 
ability of detecting at least one outburst for this sample size 
over this survey duration is exactly 90%. For higher event 
rates we would detect more outbursts and would improve our 
knowledge of the rate. 

ure[^ Access to a 100 year baseline would mean either 
a reduction in the number of stars needed in the survey, 
or an improvement on the limit able to be placed on the 
FU Ori rate - both by the same factor of 100 - relative 
to a 1 year baseline survey. 

In practice, an automated comparison of two fields at 
widely separated epochs would select stars with mea¬ 
sured magnitudes that differed by some threshold Am as 
candidate FU Ori objects, to be confirmed or rejected by 
follow-up observations. The actual magnitude increase 
AM to which such a strategy is complete is always larger 
than Am because of measurement errors in the photom¬ 
etry at each epoch. Let a star have a true magnitude Mi 
in epoch 1 and a true magnitude M 2 in epoch 2. Let the 
star have observed magnitudes mi and m 2 that can be 
modelled by Gaussian random variables with respective 
means Mi and M 2 and variances (t\ and (t| . The magni¬ 
tude difference is then also a Gaussian random variable 
with mean M 2 — My and variance + In this model, 
m 2 — mi > (M 2 — Ml) — 1.282^yaY^HT^. with 90% 
probability. A survey with a cutoff Am for identifying 
an outburst detection will be 90% complete to outbursts 
with amplitude larger than Am + 1.282-^CTj + 

In an actual data analysis, the value of Am can be cho¬ 
sen arbitrarily based on the expected outburst amplitude 
and the quality of the available photometry. Statistical 
outliers based on photometric errors are a minor con- 











taminant even if one targets only the weaker outbursts. 
For AM of 2 mag and poor photographic photometry 
such as CTi = 0.3 mag, Am = 1.46 mag with a statistical 
false positive rate of 3 x 10“"^ star“^. Raising the magni¬ 
tude threshold, for example to exclude astrophysical false 
positives drawn from the large variety of young star and 
other types variables that show shorter-term variations 
at the roughly 2 mag level, would effectively increase the 
sample size needed to ensure good statistics. 

4.1. Recent Results on High Amplitude Outbursts 

Attempting to address the FU Ori rate problem, Scholz 
(2012) used the near-infrared 2MASS-UK1DSS ~8-year 
baseline to investigate several hundred stars and Scholz 
et al. (2013) used the mid-infrared Spitzer-WISE 5- 
year baseline to compare photometry of several thousand 
stars. Consequently, as clearly illustrated in Figure 
they could constrain the FU Ori rate to only much higher 
values (by several orders of magnitude) than the histor¬ 
ical estimates of the FU Ori rate that were discussed in 
§3.1, established based on alternate considerations. 

From a survey with PTF consisting of over 1000 obser¬ 
vations over the past 5.5 years towards the North Amer¬ 
ica and Pelican Nebulae, two large-amplitude brightening 
events were discovered: PTF lOnvg (V2492 Cyg) which 
is likely a combination of accretion burst behavior and 
variable obscuration (Hillenbrand et al. 2013) perhaps 
analagous to V1647 Ori, and PTF lOqpf (V2493 Cyg) 
which is a bona fide FU Ori star (Semkov et al. 2010; 
Miller et al. 2011) previously known as the Class II clas¬ 
sical T Tauri star LkHa 188/G4or HBC 722. Both were 
discovered within the first season of PTF monitoring and 
- given the only few thousand star sample of known 
young stars in the region- would suggest an FU Ori rate 
significantly higher than the canonical value. 

Considering a parent sample of ~2000 stars in the re¬ 
gion with evidence of circumstellar material, and assum¬ 
ing perfect detection efficiency, the implied event rate is 
then 9 x 10“® outbursts per star per year. From the 
simulations presented here, the 90% confidence interval 
for this rate is [3 x 10“®,4 x 10“^] outbursts per star per 
year. Of course, statistical estimates based on only a sin¬ 
gle detection are not adequate predictors. At the implied 
event rate, the chances of not seeing another FU Ori out¬ 
burst for another 20 years are only 3%, while at the lower 
bound, there is a 30% chance of no additional outbursts 
in the next 20 years. 

4.2. Recent Results on Low Amplitude Bursts 

In addition to the bona fide FU Ori star outbursts, 
there are also lower-amplitude bursts that have been 
identihed, including the EX Lup group illustrated in Fig¬ 
ure 1, as well as even shorter timescale and lower am¬ 
plitude bursts that are more routinely found in modern 
survey data. These less extreme but more common types 
may also play an important role in the accretion history 
of a star. 

Among a sample of young stars selected based on vari¬ 
ability and infrared-excess, Findeisen et al. (2013) found 
in the North America and Pelican Nebulae region that 
12% ± 3% exhibited bursting behavior on timescales 
from 0.1 to 100 days, with amplitudes from 0.2 to 1.5 
mag. If confirmed as accretion-driven events, the cases 


of repeated bursts reported there may be the first empir¬ 
ical examples of accretion-driven cyclic behavior caused 
by processes in the inner disk, such as those predicted on 
time scales of a few tens of days in the disk instability 
mechanisms described in §2.2. 

The simulations presented above are not applicable to 
these shorter duration burst events since we have not 
included either the survey cadence (which needs to be 
short enough to catch and time-resolve the burst) or the 
burst duration (which needs to be long enough to be 
detected) as model parameters in our simple framework. 

4.3. The Numerator 

As stated earlier, estimating an event rate requires 
measuring a number of detections, with considerations 
to bias, from among the number of objects for which de¬ 
tections could be made. The numerator in the rate equa¬ 
tion is the number of FU Ori outbursts observed. There 
are in fact two sides of the FU Ori outburst phenomenon 
which may be amenable to detection: one the rise side or 
outburst itself, occuring over weeks or months, and the 
other the decay side, occuring over decades and centuries. 

Numerous projects dedicated to other (usually ex- 
traglactic) science areas are also finding young star out¬ 
bursts, some of which are even FU Ori like. The cur¬ 
rently operating ASAS, CSS, and PTF will be succeeded 
by next generation surveys such as ATLAS, ZTF, and 
LSST, all of which are capable of making discoveries in 
this area. The ATLAS project, for example, is predicting 
5 FU Ori stars in its first year of operations. 

Survey yield will depend on parameters such as sensi¬ 
tivity, cadence, and galactic latitude coverage, as well as 
on astrophysical factors such as the range in burst am¬ 
plitudes and rise times, and the frequency of visibility at 
optical vs infrared-only vs submm-only wavelengths (e.g. 
Johnstone et al 2013). None of these considerations are 
included in our simple simulations, in which they were 
wrapped up in the factor e which was set to unity. How¬ 
ever, the results may be generalized to the case of e 7^ 1 
by considering that, even if not all FU Ori sources are de¬ 
tected, a Poisson/binomial process will still be observed 
albeit at an effectively reduced rate. A best estimate rate 
and a confidence interval could still be calculated, then 
multiplied by 1/e. 

In addition to catching the FU Ori rise phase, various 
of these imminent time domain sky survey projects could 
also identify FU Ori objects during the decay phase, by 
selecting for follow-up spectroscopy sources with consis¬ 
tent fading over their 5-10 years of operations. However, 
there is significant diversity in the post-outburst light 
curves of known FU Ori stars (see e.g. Clarke et al. 
2005, their Figure 1). For example, FU Ori itself has 
decayed only slightly in 50 years, while V1057 Cyg had 
a much more rapid decay; V1515 Cyg was slower to rise, 
but has had a relatively flat post-outburst light curve. 

4.4. The Denominator: Census of Selected Star 
Forming Regions 

The demoninator in the rate equation is the number 
of young stars with the possibility to undergo an FU Ori 
outburst in a given survey. The true survey size is Nyso 
young stars with disks, which given a certain survey area 
is generally far less than as the number of objects photo¬ 
metrically detected in the field. Young stars are iden- 
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tified via a variety of mechanisms associated with ei¬ 
ther stellar or circumstellar processes. Taken in isolation, 
none may not be 100 % reliable indicators of membership 
in <1-3 Myr old star forming regions, but in combination 
they can lead to secure knowledge of the total size and 
identity of the targetted young stellar population. 

The usually considered diagnostics of young T Tauri 
stars are: photometric variability, low excitation emis¬ 
sion lines, presence of quick-burning lithium, infrared 
excess due to circumtellar dust, ultraviolet excess due 
to re-radiated accretion shocks, and x-rays from both 
accretion and stellar coronal activity. As the FU Ori 
phenomenon is driven by massive disks, only so-called 
Class-I and Class-II stars should be included in assess¬ 
ing survey statistics. This is a much easier scenario than 
other statistical problems such as the lifetime of disks 
where both the disked and non-disked populations must 
be thoroughly known, as well as the absolute source ages. 

There has been great progress over the last decade in 
terms of census building in star-forming regions. How¬ 
ever, the job is not complete - even for regions within 
only a few hundred parsecs of the Sun. Spitzer-identified 
(infrared excess) and Chandra-identified (x-ray) sam¬ 
ples are incomplete, even in combination, in part be¬ 
cause these were pointed missions. The recent availabil¬ 
ity of wide-field survey data from WISE (mid-infrared), 
2MASS and UKIDSS (near-infrared), iPHAS/UVEX 
(optical), and GALEX (ultraviolet) can help fill in these 
census gaps near and in the galactic plane. Non-outburst 
photometric variability is also - after some decades of 
hiatus - again being used to identify young stars. GAIA 
results on kinematics will help define cluster and moving 
group membership for optically visible young stars. 

5. DISCUSSION AND CONCLUSIONS 

There are at present reasonable empirical constraints 
on the rate of occurence of phenomena such as M dwarf 
flares, novae, supernovae, and gamma ray bursts. Each 
can be quantified in terms of the number per year per 
square degree per magnitude interval. The young star 
FU Ori outbursts, however, which have luminosity in¬ 
creases between those of M dwarf flares and novae, are 
rare enough that only crude estimates exist for the num¬ 
ber that occur within each young star disk lifetime. And 
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while the rate is a significant number, more meaningful 
would be knowledge of the distribution of time separa¬ 
tion between FU Ori events during the first one to few 
Myr when the star first becomes visible as a self-luminous 
hydrostatically contracting object. 

In this paper we have provided a simple framework for 
considering how well we might be able to constrain the 
rate of FU Ori outburst events. The method employs 
binomial statistics which are generally applicable to any 
phenomenon that can be treated as either occurring (e.g. 
an outburst is detected) or not occuring (no outburst is 
detected). In our presentation we have provided numbers 
based on a 90% measurement confidence requirement (a 
parameter of the statistical analysis that depends on the 
data actually obtained) with a 90% probability (a pa¬ 
rameter which refers to having a representative sampling 
of the astrophysics), under the assumption that all out¬ 
burst which occur are detected. Figures 2, 3, 4, and 5 
illustrate the simulation results. 

With the rigor specified above, empirically con¬ 
straining the FU Ori outburst rate to a level of < 
lO”"* Mq year“^ requires > 10^ young stars with disks 
to be monitored for > 1 year, or > 10^ such stars for 

> 10 years. Similarly, constraining the FU Ori out¬ 
burst rate to < 10“^ Mq year“^ requires > 10® young 
stars with disks to be monitored for > 10 years, or 

> lO"* such stars for > 100 years. Finally outburst rates 
< 10“® Mq year“^ could be probed if > 10® young stars 
with disks were monitored for > 10 years, or > 10® such 
stars for > 100 years. In any case, more than several 
lO"* suitably disked young stars are needed in order to 
be on the linear parts of these scaling relations, given 
the likely range of true outburst rates. An extension of 
our calculation to shorter duration event types such as 
EX Lup / V1647 Ori or lower level bursts would require 
consideration of the survey cadence relative to the burst 
duration as additional model parameters. 

Einally, we note that in order to estimate the true out¬ 
burst rate, attention is needed not only to the numerator 
that reflects the number of actually detected outburst 
events, but to careful definition of the denominator, that 
is, the number of comparably aged young disked stars in 
the studied region. 
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